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We study modulational instability (MI) of plane waves in nonlocal nonlinear Kerr media. For 
a focusing nonlinearity we show that, although the nonlocality tends to suppress MI, it can never 
remove it completely, irrespectively of the particular profile of the nonlocal response function. For a 
defocusing nonlinearity the stability properties depend sensitively on the response function profile: 
for a smooth profile (e.g., a Gaussian) plane waves are always stable, but MI may occur for a 
rectangular response. We also find that the reduced model for a weak nonlocality predicts MI in 
defocusing media for arbitrary response profiles, as long as the intensity exceeds a certain critical 
value. However, it appears that this regime of MI is beyond the validity of the reduced model, if 
it is to represent the weakly nonlocal limit of a general nonlocal nonlinearity, as in optics and the 
theory of Bose-Einstein condensates. 
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I. INTRODUCTION 

Modulational instability (MI) constitutes one of the 
most fundamental effects associated with wave propa- 
gation in nonlinear media. It signifies the exponential 
growth of a weak perturbation of the amplitude of the 
wave as it propagates. The gain leads to amplification of 
sidebands, which breaks up the otherwise uniform wave 
front and generates fine localized structures (filamenta- 
tion). Thus it may act as a precursor for the formation 
of bright spatial solitons. Conversely the generation of 
dark spatial solitons requires the absence of MI of the 
constant intensity background. 

The phenomena of MI has been identified and studied 
in various physical systems, such as fluids plasma [^, 
nonlinear optics discrete nonlinear systems, such as 
molecular chains 1 5 1 and Fermi-resonant interfaces and 
waveguide arrays |'], etc. It has been shown that MI is 
strongly affected by various mechanisms present in non- 
linear systems, such as higher order dispersive terms in 
the case of optical pulses [[7| , saturation of the nonlinear- 
ity and coherence properties of optical beams 

In this work we study the MI of plane waves propagat- 
ing in a nonlinear Kerr type medium with a nonlinearity 
(the refractive index change An, in nonlinear optics) that 
is a nonlocal function of the incident field. We consider 



a phenomenological model, in which the nonlocal nonlin- 
earity, induced by a wave (e.g. an optical beam) with the 
intensity I{x, z), can be represented in general form as 



An(/) 



R{x' — x)I{x , z)dx' , 



(1) 



where x is the transverse spatial coordinate and s=l 
(s=— 1) corresponds to a focusing (defocusing) nonlin- 
earity. The evolution coordinate z can be time, as for 
Bose-Einstein Condensates (BEC's), or the spatial prop- 
agation coordinate, as for optical beams. We consider 
only symmetric spatial response functions that are pos- 
itive definite and (without loss of generality) obey the 
normalization condition 



R{x)dx — 1. 



(2) 



Thus we exclude asymmetric effects, such as those gener- 
ated by asymmetric temporal response functions (with x 
being time) , as in the case of the Raman effect on optical 
pulses | |To[ |. 

In nonlinear optics Eq. (^ represents a general phe- 
nomenological model for media in which the nonlinear 
refractive index change (or polarization) induced by an 
optical beam is determined by some kind of a transport 
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process. It may include, e.g., heat conduction in mate- 
rials with a thermal nonlinearity [pl] - [l3| or diffusion of 
molecules or atoms accompanying nonlinear light prop- 
agation in atomic vapours jl4j. Nonlocality also accom- 
panies the propagation of waves in plasma [p^p^ , and a 
nonlocal response in the form (|l|) appears naturally as a 
result of many body interaction processes in the descrip- 
tion of Bose-Einstein condensates . 

The width of the response function R{x) relative to the 
width of the intensity profile I{x, z) determines the de- 
gree of nonlocality, as illustrated in Fig. |^. In the limit of 
a singular response, R(x)—6(x) [see Fig. |l|(a)], the non- 
linearity is a local function of the intensity, 

An = s/, (3) 

i.e. the refractive index change at a given point is solely 
determined by the wave amplitude at that very point. 
For this well-known Kerr nonlinearity, which appears in 
all areas of physics, MI depends only on the sign s (H. 

With increasing width of R{x) the intensity in the 
vicinity of the point x also contributes to the index 
change at that point. For a weak nonlocality, when the 
width of R(x) is finite but still small compared to the 
width of l{x,z) [see Fig. |l|(b)], one can expand I{x' ,z) 
around x'=x in (n|) and obtain the simplified model 

1 /"^ 

An = s{I + -fdll), 1^7: x^R{x)dx, (4) 

^ J -00 

where 7 is a small positive parameter. This diffusion 
type model of the nonlocal nonlinearity is a model in 
its own right in plasma physics, where 7 can take any 
sign jl^Jlj. It was also appHed to EEC's |2l[, nonhnear 
opticspj] and in the continuum limit of the theory of 
energy transfer in biomolecular systems . In contrast 
to the local Kerr limit the MI now depends not only on 
the sign s, but also on the intensity of the plane wave 
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FIG. 1. Different degrees of nonlocality, as given by the 
width of the response function R{x) and the intensity profile 
I{x). Shown is the local (a), the weakly nonlocal (b), the 
general (c), and the strongly nonlocal (d) response. 

In the limit of a strongly nonlocal response, where R{x) 
is much broader than the intensity profile [see Fig. 0(dy , 
one can expand the response function around x'—x in dl) 
and obtain the simple linear model 

An = s{cq + cix + C2x'^), (5) 

where co_2 are constants. Since this model is linear all 
plane waves are stable, - there is no MI. The evolution 
of optical beams in such a medium was considered in the 
case when ci—0 [p5| , which requires the center of mass 
of the beam to be always zero. 

Modulational instability has thus been studied in the 
different limits. However, in the general case (|^) the 
problem of MI has only been studied for a few particular 
cases of the response function R{x). Here we present an 
analytical study of the general case with arbitrary profile 
R{x) and confirm the results for three specific examples. 
We further show that the MI results of the weakly non- 
local model (^) lead to erroneous conclusions when used 
to predict the behavior of the full model (|^). 

II. GENERAL THEORY 

We will consider an optical beam propagating along 
the z-axis, with the scalar amplitude of the electric field 
expressed as 

E{x, z) = i^(x, z) exp{iKz — iQ,t) + c.c, (6) 

where K is the wavenumber, U, is the optical frequency 
and ipix^ z) is the slowly varying amplitude. Substituting 
Eqs. (|l|) and (||) into Maxwell's equations we obtain the 
nonlocal nonlinear Schrodinger (NLS) equation 

1 f°° 

id^^ + -dlijj + sijj R{x' ~ x)I{x\z)dx' ^0, (7) 

2 J -00 

where I(x, z)=\il^{x, z)|^ is the intensity of the beam. The 
model (|7|) permits plane wave solutions of the form 

ip{x, z) = ^/p^exp{ikox - iujoz), Po > 0, (8) 

where po, fco, and luq are linked through the nonlinear 
dispersion relation 

^0 = ^fco - sPo, (9) 
which is the same as for the standard local NLS equation 

id,i;+^dl^ + s\^P\''ij^O. (10) 
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Next, let us carry out a linear stability analysis of the 
plane wave solutions (H). Assume that 



a = l/(4po) 



(20) 



ill) 



where ai(a;, z) is a small complex perturbation. Inserting 
this expression into the nonlocal NLS equation (0) and 
linearizing around the solution (||) yields the evolution 
equation for the perturbation 

/•OC 

+2spo / i?(C'-ORe{ai(e',^)}'^a;' = 0. (12) 



In deriving this equation we have used the nonlinear dis- 
persion relation (^) and gone into a coordinate frame 
moving with the group velocity Cg = dcoo/dko = k^, i.e. 



T = Z, — X — CgZ 



(13) 



Decomposing the perturbation into real and imaginary 
parts, ai — u + iv, we obtain two coupled equations 



drU + \d'^v = 0, 

drv - \dlu - 2spo R{C - r)dC = 0. 

By introducing the Fourier transforms 



(14) 



J — oo 

/OC 
-oo 
/•oo 

R{k) = / i?(0exp[iA:C]d5, 



(15) 



and exploiting the convolution theorem for Fourier trans- 
forms, the linearized system is converted to a set of ordi- 
nary differential equations in fc-space 



drU — \k'^v — 0, 

drV ~\- ^k'^u — 2spoRu = 0, 
which can be written in the compact matrix form 

drX ^AX, 
where the vector X and matrix A are defined as 



X 







-^k^ + 2spoR{k) 



The eigenvalues A of the matrix A are given by 



(16) 



(17) 



(18) 



(19) 



Notice that if the response function R(x) is real and 
symmetric then so is the Fourier spectrum of R{x), i.e. 
Rik) =R*ik)=R{-k). 

The general eigenvalue equation (|l^) constitutes the 
basis of our further study. First of all we notice that 
i?(0)=l, since the response function is normalized ac- 
cording to Eq. (H). For a focusing nonlinearity (s=-|-l) we 
therefore always have, by continuity, that ak^ — sR{k) < 
0, and thus A^ > 0, in a certain wavenumber band sym- 
metrically centered about the origin, where k is suffi- 
ciently small. It follows that we always will have (long 
wave) MI in the focusing case, independently of the de- 
tails in the behavior of the response function. 

In contrast, in the defocusing case (s=— 1), the stabil- 
ity properties depend in a sensitive way on the response 
function. Thus the nonlocality could possibly lead to MI, 
which would never be present in the defocusing local NLS 
equation (p^. 

The well-known modulational instability (stability) re- 
sult for the standard local NLS equation (|o|) is easily 
recovered from the general eigenvalue equation (^9|) by 
setting R{x)=6{x), where S{x) is the Dirac delta func- 
tion: We get 

X"^ = -k'^poiak'^ - s) (21) 
where s = +1 (s = — 1) yields instability (stability). 



III. WEAKLY NONLOCAL LIMIT 

The MI properties in the local limit (^), described by 
the NLS equation (|o|), are well-known and the strongly 
nonlocal limit (^ is linear, so obviously it does not dis- 
play MI. Here we briefly consider the interesting regime 
of a weak nonlocality, i.e., the situation when the typical 
width of the response kernel is small compared to the 
characteristic length of the wave modulation (see Fig. |^). 
In this case the nonlinearity has the form (H) and the 
NLS equation (|) reads 



(22) 



where we have defined the parameter a as 



where the nonlocality or diffusion parameter 7 is defined 
in Eq. Equation ( |2^ ) has been discussed in the lit- 
erature in the context of plasma physics [0 and the 
continuum limit of discrete molecular structures ^2 23 
and has been shown to posses bright and dark soliton 
solutions. Their exact analytical soliton form was found 
recently [Q. Importantly, 7 is small and positive when 
Eq. ( p^ ) represents the weakly nonlocal limit of the gen- 
eral nonlocal model (|^), whereas 7 can take any sign in 
plasma physics, where Eq. ( |2^ ) is a model in its own 
right. 
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The instability condition for the weak nonlocality is 
most easily obtained by Taylor expanding the spectrum 
of the response kernel about k=0 and only retain the two 
lowest order nontrivial terms in the expansion, i.e. 

^(A;)~l-7fc2, (23) 

where obviously 7fc^ ^ 1 must be fulfilled. The eigen- 
value equation then reads 

A2 = -fc2^o[(a + s7)fc2_s]. (24) 

This result does not depend on the detailed behavior of 
the response kernel, as should be expected. 

The eigenvalue Eq. (^ ) clearly shows that MI is al- 
ways present for a focusing nonlincarity (s=+l), since 
for sufficiently small wavenumbers > for all values 
of a and 7. This is in correspondance with the results of 
the general nonlocal model (0). 

More interesting is the defocusing case (s=— 1), where 
there is a critical value of the intensity, 

Po = 1/(47), (25) 

which depends on the degree of nonlocality through the 
parameter 7, and corresponds to a = 7. When 7 is posi- 
tive (and s = —1) low intensity plane waves with pQ < p'^ 
are modulationally stable, whereas MI appears at high 
intensities, po > , despite the nonlinearity being de- 
focusing. This MI is short-wave, appearing for pertuba- 
tions with wavenumbers exceeding the threshold value 



k" = (26) 

The stability criterion po < Po' corresponds exactly to 
the criterion for existence of dark soliton solutions in 
this case [Q. The existence of MI in the weakly non- 
local model is in sharp contrast to the local defocusing 
Kerr nonlinearity, for which there is no MI. 

In the case 7 < in Eq. (p2[) , relevant in plasma 
physics, nonlocality does not affect the MI properties: 
For a focusing (defocusing) nonlinearity all plane waves 
are modulationally unstable (stable). In Table. I we 
have summarized the MI results of the weak nonlocal- 
ity model. 
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TABLE I. Modulational stability properties of the 
weakly nonlocal model (^3). 
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FIG. 2. Growth rate vs. wavenumber for plane waves with 
intensity po = l in media with a weakly nonlocal focusing (s=l, 
solid lines) and defocusing (s=-l, dashed lines) nonlinearity. 

So only for defocusing media with 7 > does a weak 
nonlocality change the MI properties in terms of exis- 
tence, as compared to the local NLS limit. However, 
nonlocality does always affect the structure of the MI 
gain spectrum, as illustrated in Fig. || for 7 > 0. In the 
focusing case Fig. ^ clearly shows how the nonlocality 
tends to suppress MI by decreasing both the maximum 
gain and the gain bandwidth. In contrast, in the defo- 
cusing case, the nonlocality promotes MI, - the higher the 
degree of nonlocality 7 the lower the threshold value of 
the wavenumber, above which the instability develops. 




FIG. 3. Development of MI of a plane wave with intensity 
po = l in a weakly nonlocal defocusing (s=-l) medium. The 
initial amplitude is modulated with the wave vector k=7Tv. 

In Fig. 3 we show the development of MI of a plane 
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wave in a self-defocusing weakly nonlocal medium, as 
obtained by numerical integration of the NLS equation 
(p^). As an initial condition we used a plane wave with 
periodically perturbed amplitude. It is evident that this 
perturbation grows as the wave propagates. 

Although Eqs.dl) and (|6|) predict MI in the defocus- 
ing case it is important to stress that this result does 
not actually apply to the case of weak nonlocality. In 
deriving the model (|2^) and the expansion ( p3| ) it was 
assumed that 



-yP < 1, 



(27) 



which means that the width of the response function 
must be small compared to the characteristic length-scale 
of the modulation. However, from Eq. ( p^ ) we get that 
MI occurs when 



{a - -i)k^ + 1 < 0, 
which can be rewritten as 



7fc^ > 1 + ak^. 



(28) 



(29) 



Since a — l/(4po) is positive this condition cannot be 
fulfilled without violating the perturbation requirement 
( p7| ) . Hence, in the region of parameters where MI occurs 
(and dark solitons do not exist), the weak nonlocality ap- 
proximation breaks down. This is in accordance with the 
analysis of the full nonlocal model (0), which, as we show 
below, does not display MI for self-defocusing nonlinear- 
ities with a smooth nonlocal response function. 

It should be stressed however that the appearance of 
MI is a valid result in the situation when the weakly non- 
local model (^ is a model in its own right, as for example 
for upper hybrid modes in plasma p7|]. 



IV. FULLY NONLOCAL RESPONSE - 
EXAMPLES 

We will now discuss the stability issue in the general 
nonlocal Kerr nonlinearity with an arbitrary degree of 
nonlocality by considering a few specific examples of the 
nonlocal response function. 



A. Gaussian response function 

Let the response function of the nonlocal nonlinear 
medium have the following Gaussian form 



R{x) 



1 



exp 



„2 1 



(30) 



In this case 7 = cr^/2, the Fourier transform of the re- 
sponse function is given by 



R{k) ~ exp 



1 



(31) 



and the eigenvalue equation ( pj| ) reduces to 



= —k^po I ak'^ — scxp 



— a k 
4 



(32) 



For the defocusing case (s = — 1) we find modulational 
stability for all values of the intensity po, in contrast to 
the behaviour predicted by the weakly nonlocal model 
(p^, according to which plane waves become unstable 
if their intensity is higher than the critical intensity p". 
Thus, as expected, this proves that the weakly nonlocal 
model gives incorrect predictions of the MI properties of 
the full general nonlocal model. 

A simple analysis reveals finite band width MI in the 
focusing case (s = -1-1) with a growth rate given by 



|ReA| = ^pJ^Jexp 



(33) 



The dependence of the growth rate on the wave num- 
ber and degree of nonlocality a is displayed in Fig. ^. 
Evidently, as the degree of the nonlocality increases the 
MI bandwidth shrinks and the maximum growth rate de- 
creases. Thus nonlocality clearly tends to suppress MI in 
this also predicted by the weakly nonlocal model. 
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FIG. 4. MI gain profiles for po=l in self-focusing {s- 
nonlocal Kerr media with a Gaussian response profile. 
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In Fig. H we demonstrate the influence of nonlocality 
on the dynamic development of MI. To this end we inte- 
grated numerically the NLS equation with the Gaus- 
sian response function (^). As an initial conditions we 
used a unit intensity plane wave with imposed weak pe- 
riodic perturbation 



4>{x, z^O)^l + IQ-'^ cos(1.5a;) 



(34) 



It is evident that increasing the degree of nonlocality a 
suppresses the instability by lowering the gain. 




FIG. 5. Development of MI of a unit amplitude plane wave 
in a nonlocal self-focusing (s=l)medium with a Gaussian re- 
sponse function for (a) a=0.1 and (b) a=1.0 



B. Exponential response function 

Another interesting and important example involves 
the exponential response function 



R{x) 



1 

2^ 



exp 



(35) 



for which 7 = 2a^ , with a again determining the degree of 
nonlocality. The importance of this particular response 
lies in the fact that it applies exactly to media with a non- 
linearity (i.e. nonlinear refractive index An{x)) described 
by the diffusion equation 



An = IV'I^ 



(36) 



The equivalence between the nonlocal nonlinearity (|i|) 
with the exponential response ( |35| ) and the diffusion 
model (^6|) is easily shown by Fourier transformation. 

As discussed by Litvak et al. jl^ the model (^ ) de- 
scribes the nonlinear thermal self-focusing of electromag- 
netic waves in a weakly ionized plasma. In this particular 
case An corresponds to the change of the plasma tem- 
perature. From the Fourier transform 



R{k) 



1 



1 



we obtain the eigenvalue equation 



1 + (crfc)2 



(37) 



(38) 



which coincides with that derived in Ref. flsl] [Eq. (4.1) 
in this reference]. The dispersion relation (|38|) imposes 
MI properties similar to those found for the Gaussian 
response. In particular, it follows that in the self- 
defocusing regime (s = —1) the plane wave solution is 
always modulationally stable. On the other hand, in the 
self-focusing case plane wave develops transverse insta- 
bility with the growth rate given by 



|ReA| = \k\^. 



1 



{aky 



(39) 



In Fig. (^ we show the MI growth rate versus the wave 
number for a few values of the nonlocality parameter a. 
As for the Gaussian response the nonlocality tends to 
suppress MI by narrowing the gain bandwidth and de- 
creasing the maximum gain. However, only in the ex- 
treme limit of an infinitely nonlocal response (cr — > 00) 
does MI completely disappear, as the medium no longer 
exhibits any nonlinearity. 
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FIG. 6. MI gain profiles for po=l in nonlocal self-focusing 
(s=l) media with an exponential response function. 
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C. Rectangular response function 

Finally we consider the rectangular response function 

r 5L for |x| < a 
R{x) = <^ (40) 



for \x\ > (J 



for which 7 = cr^/3. From the Fourier transform 



R(k) = 



sm{ka) 
ka 



we obtain the eigenvalue equation 

sin(fc(T) 



ak^ 



ka 



(41) 



(42) 



In the focusing case one again finds that there is always 
(for any value of a and a) a finite bandwidth MI gain 
band symmetrically located about the origin in /c-space, 
i.e. for I A: I < fccr, where 



2 _ sin ( fccr g) 



(43) 



This is illustrated in Fig. 7, where we show the gain pro- 
files for several degrees of nonlocality a. Again MI is 
being suppressed by the nonlocality, but now additional 
MI bands appear for large values of cr (see (t—20). 
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FIG. 7. MI gain profiles for po~l in self-focusing (s=l) 
nonlocal media with a rectangular response function. 

In the defocusing case a careful analysis reveals that 
for large and moderate values of the parameter a/tr^, 
i.e. small and moderate values of poc'^, plane waves are 
modulationally stable, while in the opposite regime one 
or more MI bands appear, as shown in Fig. 8. This insta- 
bility is caused by the negative-sign bands of the sinc(a;)- 
profile Fourier Transform (^l|), with the first instability 



band being in the region n < ka < 27r. More specifically 
the instability appears when the background intensity 
exceeds the threshold value 



p*o > 21.05/cr^ 



(44) 
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FIG. 8. MI gain profiles for po=l in self-defocusing (s=-l) 
nonlocal media with a rectangular response function. 



V. CONCLUSIONS 

We studied the modulational stability properties of 
plane waves in nonlinear media with a general nonlocal 
Kerr nonlinearity, described by a nonlocal NLS equation. 
We derived an analytical expression for the growth rate 
of the instability for arbitrary response profiles, i.e. for 
arbitrary degrees of nonlocality. We found that in the 
self-focusing case the plane wave solution is always un- 
stable, - the nonlocality tends to suppress the instability, 
independently of the form of the response function, but 
can never suppress is completely. 

In self-defocusing media the stability properties de- 
pend strongly on the specific form of the nonlocal re- 
sponse. For any strictly positive continuous response 
profile all plane waves are stable for any degree of nonlo- 
cality. However, exotic profiles with jump-discontinuities 
like the rectangular response, may actually promote MI 
when the intensity is sufRciently high" 

In the so-called weakly nonlocal limit (narrow response 
function compared to the intensity profile) we find that 
plane waves become modulationally unstable if their in- 
tensity exceeds a certain critical value. However, we show 
that this surprising prediction of MI in the defocusing 
case can never apply to the full nonlocal model with a 
smooth response profile. Thus one should be cautious 
with applying results of the simple weakly nonlocal model 
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to the real general case, in particular when studying MI 
and dark solitons. 
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